Linear and Nonlinear PT-symmetric Oligomers: A Dynamical Systems Analysis 
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In the present work we focus on the cases of two-site (dimer) and three-site (trimer) configura- 
tions, i.e. oligomers, respecting the parity-time (PT) symmetry, i.e., with a spatially odd gain-loss 
profile. We examine different types of solutions of such configurations with linear and nonlinear 
gain/loss profiles. Solutions beyond the linear PT-symmetry critical point as well as solutions with 
asymmetric linearization eigenvalues are found in both the nonlinear dimer and trimer. The latter 
feature is absent in linear PT-symmetric trimers, while both of them are absent in linear PT- 
symmetric dimers. Furthermore, nonlinear gain/loss terms enable the existence of both symmetric 
and asymmetric solution profiles (and of bifurcations between them), while only symmetric solutions 
are present in the linear PT-symmetric dimers and trimers. The linear stability analysis around 
the obtained solutions is discussed and their dynamical evolution is explored by means of direct 
numerical simulations. Finally, a brief discussion is also given of recent progress in the context of 
PT-symmetric quadrimers. 



I. INTRODUCTION 



In the late 1990s, a radical yet well physically motivated proposal emerged in the context of the study of funda- 
mentals of quantum mechanics. This was the suggestion of Bender and co-workers [l[ that Hamiltonians that respect 
(— i . two principal physical symmetries of the dynamics, namely Parity (P) and Time-reversal (T) could enable the iden- 
tification of real eigenvalues, which is a property which is highly desirable for operators associated with measurable 
quantities, even if the associated Hamiltonians are not Hermitian. A caveat in that regard, however, was that as 
the gain/loss parameter introduced in these non-Hermitian Hamiltonians was varied, a transition was quantified and 
termed the PT-symmetry breaking transition, beyond which the eigenvalues of the relevant operator were no longer 
purely real. It is interesting to note here that a prototypical playground where such ideas can be explored is that of 
standard Schrodinger Hamiltonians of the form H = — (1/2)A + V(x), for which in the case of a complex potential 
V, it is straightforward to see that the above constraints of PT symmetry amount to the potential satisfying the 
£^ ; condition V(x) = V*(-x). 

While the initial proposal of the above possibility was one of mathematical origin (rather than one inspired by a 
specific physical application), in the past few years, a significant vein of potential applications of such Hamiltonians has 
been initiated, predominantly so in the field of nonlinear optics. There, the work of Christodoulides and co-workers Q 
gave rise to the realization that the synthetic systems that can be engineered therein enable a potential balance of 
gain (through suitable amplifiers) and the abundantly present losses in order to produce experimental realizations 
of PT-symmetric systems. An additional feature present in such settings which made both their theoretical and 
experimental investigation even more interesting was the presence of nonlinearity which, in turn, rendered worthwhile 
the exploration of the dynamics of nonlinear waves (such as bright or gap solitons [|| and more recently of dark solitons 
and vortices (3[). The above optical settings were in fact the ones that enabled the first experimental realizations 
of PT-symmetry. This was done in the context of waveguide couplers (i.e., either two waveguides with and without 
loss [4| - the so-called passive PT- or in the more "standard" case of one waveguide with gain and one with loss [5j)- 
More recently, electronic analogs of such systems have been engineered in the work of [f|. In parallel to these 
pioneering steps in the realm of experiments, numerous theoretical investigations have arisen both in the context 
of linear PT-symmetric potentials [7H23j and even in the case of the so-called nonlinear PT-symmetric potentials 
(whereby a PT-symmetric type of gain/loss pattern appears in the nonlinear term) [2^ - [26| . 

In the present work, we revisit this theme of both linear and nonlinear PT-symmetric potentials. We do so in the 
special context of the recently proposed PT-symmetric "oligomers" , namely few site configurations, which are not 
necessarily couplers (i.e., dimers) but rather can also be trimers, quadrimers etc. For the sake of illustration, we 
restrict ourselves to dimers and trimers herein but briefly also touch upon recent works in the context of quadrimers. 
In what follows below, we study the case of nonlinear-PT-symmetric dimers [24j and trimers and subsequently restrict 
considerations to the case of their linear-PT-symmetric analogs [l3[ . One of the fundamental side-effects of the fact 
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that the nonlinearity does not commute with the PT operator is the existence of nonlinear solutions that persist 
past the linear PT-phase-transition threshold. Furthermore, one of the principal consequences of the presence of 
gain/loss terms in the nonlinearity is the existence of both symmetric and asymmetric (in their amplitude) stationary 
solutions, with the latter possessing a non-symmetric linearization spectrum. Interesting bifurcation phenomena (such 
as spontaneous symmetry breakings) are, additionally, found to arise in this case. Our presentation is structured as 
follows: we first examine the nonlinear-PT-symmetric dimer in section II (existence of solutions in II. A, linear stability 
setup in II. B and numerical results in II. C), while the corresponding analysis for the trimer is done in section III. In 
section IV, we review the corresponding case of linear-PT-symmetric oligomers, while in section V, we summarize our 
findings and present our conclusions. 

II. ANALYSIS OF STATIONARY SOLUTIONS FOR THE NONLINEAR-PT-SYMMETRIC DIMER 

CASE 

We first consider the so-called PT-symmetric coupler or dimer, in which a gain/loss pattern appears both in the 
linear and nonlinear terms. The dynamical equations have the form: 

iu t = —ev + (p r - ip lm )\u\ 2 u + iju 

iv t — -eu + (p r + ip lm )\v\ 2 v - ijv. (1) 

The model contains the Kerr nonlinearity which is relevant to optical waveguides and is effectively a generalization 
of the experimental framework of [5[, in that nonlinear (i.e., amplitude-dependent) gain and loss processes are taken 
into account. We have used u(t) and v(t) to denote the two complex- valued variables for the dimer and the evolution 
variable is t (in optics, this is actually a spatial variable standing for the propagation distance along the optical 
crystal). Considering the prototypical stationary solutions of the system, we let u(t) and v(t) have the forms: 

u{t) = ae lEt , v(t)=be lEt (2) 

where E is the propagation constant while the complex numbers a and b denote the amplitudes of the dimer sites. 
Plugging this ansatz into Eq. (fTJ), one finds the complex nonlinear algebraic equations: 

Ea — eb — (p r — ipi m )\a\ 2 a — i^a 

Eb = ea - (p r + ip im )\b\ 2 b + ijb. (3) 
We now use a polar decomposition of a and b of the form: 

a = Ae 1 ^, b = Be^ b (4) 

for real- valued A, B, <p a and 4>b- Plugging Eq. (j4]) into Eq. ([3]) and writing these equations in terms of their real and 
imaginary parts, we find: 

EA = eBcos(0 b - 4> a ) - p r A 3 
EB = eAcos((f> b - <fi a ) - p r B 3 
-eA sin(</> 6 - <f> a ) - p m B 3 +jB = (5) 
eB sin(0f, - <f> a ) + p lm A 3 - jA = 0. 

The last two equations yield 

(A 2 - B 2 ) [p. m (A 2 + B 2 ) - 7 ] =0. (6) 

We note that Eq. ([B]) yields a simple algebraic condition which connects the amplitude of the two dimer sites. This 
allows us to distinguish several subcases of interest. We look for nontrivial solutions A and B in each of the subcases 
presented in the following section. 
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A. Existence of Localized Modes for the Dimer Case 

Eq. (|6]) identifies the different scenarios for the values of A and B. We now examine the three cases that arise from 
this equation for our dimer dynamical system. 

• Case I : A 2 — B 2 and A 2 + B 2 ^ j/p im : 
Recall the equations given in ([5]): 

EA = eBcos(cj) b - (p a ) ~ p r A 3 ; EB — eAcos((p b — cp a ) — p r B 3 (7) 

and 

eB sin(0 fc -</>„) + PimA 3 - j A = 0; - eA sin(0 b - (f> a ) - p m B 3 +-/B = 0. (8) 

Since A — B (i.e., these are symmetric solutions) in this case, the two equations in each set are equivalent. 
Thus, we have: 

■(a. a\ -(.P«nA 2 --f) Pr A 2 +E 

sm(0fe - 4> a ) = , COs(</> 6 - <p a ) = . (9) 

e e 

We use the relation sin 2 (^>f, — <p a ) + cos 2 ((pi, — (p a ) = 1 to determine the following quadratic equation for A 2 : 

(p 2 r + p 2 m )A i + 2(Ep r - 1Pim )A 2 + 7 2 + E 2 - e 2 = 0. (10) 
The solution of the resulting bi-quadratic equation reads: 



B 2 = ~{Epr - IPim) ± V(E P r - 1 Plm ? = (pj + P j m )(l 2 + g = 6 2 ) ( . 

P 2 r+PL 



with the restriction that 



(Ep r - 1Pm y > (p 2 r + p 2 m )(Y + E 2 - e 2 ). (12) 



Case II : A 2 + B 2 = 7 /p tm and A 2 ^ B 
Under these conditions, one can get 



with the restriction that 



7 , / 7 2 



A = tt—±\it* 5 — ( 13 ) 



j2 7 _ / 7 2 



B = v ^i^^TT^ (14) 



E = -™ (15) 



cos((p b -(p a ) = £Pr , (16) 



VpJ+p 



.2 



7 2 e 2 

tV > -2 — r- (17) 

Pirn Pr ' Pim 



A fundamental difference of this case from case I is that here E is no longer a free parameter [24| . The solutions 
with the different amplitudes will be called asymmetric in what follows. 

• Case III : A 2 + B 2 = >y/p im and A 2 = B 2 : 
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As a final "mixed" possibility, between the above symmetric and asymmetric cases, from Eq. ((5|), it is straight- 
forward to obtain 



A = B 



7 



(18) 



with the restriction that 



COS(0 b - (p a ) 
sm((j) b - (j> a ) 

2p im E + 7/y 



2p im E + -fp r 



2epi 



J_ 

2e' 



2e 



= 1. 



(19) 
(20) 

(21) 



Once again this implies that once other parameters (such as 7, Pi m , p r and e are determined, E is not a free parameter 
but rather is obtained from Eq. (|2"TT) . These will be referred to as special symmetric solutions in the following. 

It is particularly important to highlight that both solutions of Case II (asymmetric) and ones of Case III (special 
symmetric) are present due to competing effects of the linear and nonlinear gain loss profiles; notice the opposite 
signs thereof in Eq. ((T|) and the necessity of jpi m > for such solutions to exist. In the case, where the linear and 
nonlinear gain/loss cooperate (rather than compete) such solutions would obviously be absent and the system would 
be inherently less wealthy in its potential dynamics. This point was also discussed in [24] . 



B. Linear Stability Analysis for the Dimer Case 



We now go back to our original PT-symmetric dimer with linear and nonlinear gain and loss in Eq. (JT|) and examine 
the linear stability of the solutions to this equation. We begin by setting 

u(t) = e lEt [a + pe xt + Pe yt ], v(t) = e lEt [b + qe xt + Qe A **] (22) 

where A is a complex- valued eigenvalue parameter revealing the growth (instability) or oscillation (stability) of all the 
modes of linearization of the dimer system; * denotes the complex conjugate and p, P, q, Q are perturbations to the 
solutions of interest. Plugging Eq. (|22p into Eq. (JTJ and taking only the linear terms in p, P, q and Q, we find the 
following eigenvalue problem: 

AX = iAX (23) 
where X = (p, — P*, q, — Q*) T and A is written as: 



an -a 2 {p r -ipim) -e 
{a*) 2 (p r + ip lm ) a 2 2 

-e a 33 

e (b*) 2 {p r - ipim) 



\ 



-b 2 (p r + ipim) 



(24) 



044 



/ 



where 



an = E + 2\a\ 2 {p r - ip im ) + «7 



a 22 = -2\a\ 2 {p r + ip vm ) - E + i-f 
a33 = E - 17 + 2\b\ 2 (p r + ip im ) 



a 4 , 4 , = -E-i'y-2\b\ 2 (p r -ipi m ). (25) 

The use of the symmetric, asymmetric or mixed solutions of the previous subsection into these matrix elements 
produces a 4 x 4 complex matrix whose eigenvalues will determine the spectral stability of the corresponding nonlinear 
solution. The existence of eigenvalues with positive real part X r > amounts to a dynamical instability of the relevant 
solution, while in the case where all 4 eigenvalues have A r < 0, the solution is linearly stable. 
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C. Numerical Results for the Dimer Case 

Fig. [T] shows the profile of the different branches for the dimer case and for parameters e = 1, E = 1, p T = — 2 
and p im = 1 (unless noted otherwise). The branches denoted by blue stars and red diamonds correspond to the 
case I of symmetric solutions; these two branches collide and disappear at the critical point 7 = 1.61 (when Eq. (fT2")) 
becomes an equality). The green circle and magenta cross branches correspond to case II; and the black squares 
branch corresponds to case III. For the latter two branches, when 7 is varied, E is also varied too (rather than staying 
fixed at E = 1 as for case I) according to Eqs. (fT5j) and ([2"T]) . respectively. Similar notations are used in Fig. [2j 
which shows the linear stability eigenvalues A = A r + i\ of the linearization. While the branches of case I are stable, 
it is interesting to note that the branch of case III (black squares) is stable until a pitchfork (symmetry breaking) 
bifurcation arises at 7 = 0.895 (when Eq. (p~7j) becomes an equality) and acquires a real pair of eigenvalues thereafter 
signalling its dynamical instability. On the other hand, it is at that critical point that the two branches belonging 
to case II arise. While the special symmetric black squares' branch of case III persists up to the critical point of 
7 = 2e = 2 of Eq. (|20)) . it should be pointed out that nonlinearity enables the asymmetric branches of Case II to 
persist for large values of 7, in fact well past the point of the linear PT phase transition. This feature has been 
highlighted in a number of recent works [T3, H?} ; in the case of the dimer the linear critical PT phase transition point 
is identified as 7 = e, while for the trimer setting considered below it is 7 = y/2e. 

An additional point worthy to mention here is that in linear PT-symmetric chains (just as is the case in typical 
Hamiltonian systems), if A is an eigenvalue to the linearization problem around a solution, so are —A, A, and —A 
(where the overbar denotes complex conjugation here). However, in our nonlinear PT-symmetric dimer —A and —A 
may not appear in the linearization around a particular branch, as is observed in Fig. [2] Eigenvalues of the green 
circles and magenta cross branches are not symmetric about the imaginary axis, but are symmetric with respect to 
each other. One can see from Fig. [3] that the green circles branch is always stable, while the magenta crosses branch 
is always unstable (due to an oscillatory instability associated with a complex eigenvalue pair). This is because the 
existence of asymmetry in these solutions of case II creates, in turn, asymmetries in the linearization matrix, due to the 
nonlinear gain/loss term, which breaks the PT symmetry of the linearization matrix and produces the corresponding 
observable asymmetry in eigenvalues. 

The dynamical evolution of the different elements of the bifurcation diagram of nonlinear-PT-symmetric dimer is 
shown in Fig.[3]at a fixed 7 = 1.5. In all the cases here and below, where a stationary solution exists for the parameter 
value for which it is initialized, a numerically exact solution up to 10~ 8 is typically used as an initial condition in the 
system. The system is sufficiently sensitive to dynamical instabilities that even the amplification of roundoff errors are 
enough to observe them. The stability of the case I branches is evident in the invariance of the relevant states during 
the course of the simulation (blue stars and red diamonds). On the other hand, the black squares branch is attracted 
towards the asymmetric (yet stable, as is evident in the corresponding simulation) green circles branch. Finally, the 
asymmetric magenta crosses branch leads to indefinite growth of the site with the larger amplitude (nonlinear gain) 
and the decay of the site with the smaller amplitude (nonlinear loss). 

III. ANALYSIS OF STATIONARY SOLUTIONS FOR THE NONLINEAR-PT-SYMMETRIC TRIMER 

CASE 

We now consider the generalization of the above considerations to the case of a so-called PT-symmetric trimer. 
Here, the dynamical system associated with a potential application of a three-waveguide setting is of the form: 

iu t = -ev + (p r - ip lm )\u\ 2 u + iju 
ivt — —e(u + w) — \v \ 2 v 
iw t = —ev + (p r + ip im )\w\ 2 w - i"fw (26) 



Such configurations have been considered earlier in optical applications theoretically 28j and even experimentally [29j 
in the absence of gain/loss. Here in the spirit, of @ (and also of [Hj]), we examine this case with both linear and 
nonlinear gain/loss profiles. Once again, as in the case of the dimer, we present the richer phenomenology setting 
of direct competition between linear and nonlinear gain/loss. The middle site is assumed as devoid of gain and 
loss. The Kerr nonlinearity is also assumed to be present in all three sites. Here, we use u(t), v(t) and w(t) as the 
complex-valued components for the trimer. For the stationary solutions, we again assume: 

u{t) = ae lE \ v{t)=be lEt and w(t) = ce lEt . (27) 




FIG. 1: The solution profiles of the nonlinear PT-symmetric dimer case with e = 1, p r = —2 and pi m — 1. The four panels here 
present the continuation of each branch (the amplitudes in the top left, the phases in the top right, and the real -bottom left- 
and imaginary -bottom right- parts of the linear stability eigenvalues) starting from the conservative system at 7 = 0. The five 
branches are denoted by curves of blue stars, red diamonds, black squares, green circles and magenta crosses. Blue stars: Case I 
with "-" in the amplitude; Red diamonds: Case I with "+" in the amplitude; Green circles: Case II with "+" in the amplitude 
(of A); Magenta crosses: Case II with "-" in the amplitude (of A); Black squares: Case III. Notice that the eigenvalues of green 
circles and magenta crosses are opposite to each other (see the relevant discussion in the text). We always set E — 1 in the 
case I branches, namely the blue stars and the red diamonds, which terminate at the same point when 7 = 1.61. The black 
squares are subject to a destabilizing supercritical pitchfork bifurcation at 7 = 0.895, E = 1.789 whereby the green circles and 
magenta crosses arise. The black squares branch terminates at 7 = 2; the green circles and magenta crosses exist for arbitrary 
values of the (linear) gain/loss past the linear PT-symmetry breaking point. 
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FIG. 2: The eigenvalue plots illustrating the linear stability of the nonlinear-PT-symmetric dimer with e = 1, p r = —2 and 
Pim ~ 1. For the blues stars and red diamonds branches, we use E = 1 here, while for the case II (green circles and magenta 
crosses) and case III (black squares), E is determined from the remaining parameters based on Eqs. (|15p and (|21[1 . respectively. 

Plugging Eqs. ([27} into Eq. ([55]). we find: 

Ea = eb — (p r — ipim)\a\ 2 (i — ijo,, Eb = e(a + c) + \b\ 2 b 

Ec = eb - (p r + ipim)\c\ 2 c + vyc. (28) 
Since a, b and c are complex-valued functions, we use the polar decomposition: 

a = Ae l<t, % b = Be l4,b and c^Ce 1 ^ (29) 



7 



0.51 



0.5 



0.49i 



20 
t 



40 



0.91 



0.9 



0.89 







10 
t 



20 




(a) blue stars branch 

1.5 



0.5 



(b) red diamonds branch 

3 



10 20 

t 

(c) black squares branch 



CvJ 



50 
t 

(d) green circles branch 



100 




(e) magenta crosses branch 



FIG. 3: The dynamical evolution plots of the branches for the case of the nonlinear-PT-symmetric dimer with the same 
parameter settings as in Fig. [2] when 7 = 1.5. The symmetric blue stars and red diamonds of Case I and the asymmetric green 
circles of Case II are stable, while the black squares of Case III (past the pitchfork point) and magenta crosses of Case II are 
unstable and deviate from their initial profile during the dynamics (see also the discussion in the text). 

where A,B,C,4> a ,(j) b and C are real- valued. Plugging Eq. (I2U1) into Eq. ([251) and separating the real and imaginary 
parts, we derive the following set of real- valued equations for A,B and C: 

EA = e.Bcos{(j> b - tj> a ) - p r A 3 , eBsin(0 fc - tj> a ) + p im A 3 - jA = 



EB = eAcos(cf> b — a ) + eC cos(0{, — C ) + B 3 , — eAsm(4> b — 4> a ) — eCsin( 



-.) = 



EC = eB cos{(/) b - 4> c ) - p r C 3 , eB sin(0 b - C ) - p lm C 3 + 7C = 0. 
We seek nontrivial solutions to Eqs. (15(H) . i..e, (A,B,C) 7^ (0,0,0). We can reduce Eqs. (|30l) to the form: 

( . , s EA + Pr A 3 EC + Pr C 3 

COS(0 6 -0 a ) = , COS(0 h - <p c ) = ■ 



sin((/>6 - 4> a ) 



eB 

7 A - p lm A 3 
eB 



sin (06 



:) = - 



cB 

( 7 C- Pmi C 3 ) 
eB 



Acos((j) b - (f) a ) + Ccos((j> b - 0c) = , Asin(0 fc - <j> a ) + Csin(0 b - C ) = 0. 



We apply the first four equations of Eqs. (|3Tj) into the last two equations and obtain the following relations: 
B 4 - EB 2 + E(A 2 + C 2 ) + Pr (A 4 + C 4 ) = 0, (A 2 - C 2 )[ 7 - Pm (A 2 + C 2 )] = 0. 



(30) 



(31) 



(32) 



We now determine A, B and C for several subcases (symmetric, asymmetric and mixed) as was done for the dimer 
case in section II. 
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A. Existence of Localized Modes for the Trimer Case 

For the trimer case, the special cases that can be seen to emerge for the solutions of Eqs. (|32|) can be classified as 
follows: 

• Case I : A 2 = C 2 and A 2 + C 2 _£ ~f/p im : 
In this case the algebraic equations assume the form: 

(A A ^ EA + . , 

COS(0 6 - (pa) = 5 = COS(0 fc - 4> c ) 

€13 

sin((p b -(p a ) = — — P " nA = -sin((p b - (p c ). (33) 
eB 

We now use Eq. ff32]) and cos 2 (</>{, — (p a ) + sin 2 ((pb — 4> a ) = 1 to determine: 

{Pr + pL) aS + 2 ( E Pr - lP«n)A i + (E 2 + ^ 2 )A 2 - e 2 B 2 = 0, B 4 - EB 2 + 2EA 2 + 2p r A 4 = 0. (34) 

One can solve eqns. (13"4"|) for A 2 and B 2 to complete the calculation of the relevant symmetric branch of solutions 
of Case I. 

. Case II : A 2 + C 2 = j/p im and A 2 ^ C 2 : 



From Eq. (|31[) . we obtain the four algebraic equations: 

A 2 (E + Pl A 2 ) 2 + A 2 ( 7 - p lm A 2 ) 2 = e 2 B 2 



C 2 (E + Pl C 2 ) 2 +C 2 { 1 - Pim C 2 ) 2 = e 2 B 2 

A 2 + C 2 = ^- 

Pim 

B 4 — EB 2 + E (A 2 + C 2 ) + Pr (A 4 + C 4 ) = 0. (35) 

We now have four equations but with only three unknowns (A, B and C). Therefore, in contrast to the previous 
symmetric branch of case I, one of the parameters E, e, p ri pi m , 7 is determined by the other four; i.e., not 
all four of these parameters can be picked independently in order to give rise to a solution of the trimer. Once 
again, we should nevertheless, highlight here that these asymmetric solutions only exist because of the interplay 
of linear gain/loss and nonlinear loss/gain profiles. 



Case III : A 2 + C 2 = j/p im and A 2 = C 2 : 



In this mixed case, we have 



A 2 = C 2 = 7 



with the restriction that 



B 4 - EB 2 + 2EA 2 + 2p r A 4 = 0, (36) 



£ 2 -^-^!>0. (37) 

Pim Pi m 



One can solve the following equations for B and E 



B 4 - EB 2 +2E— — + 2^-^- = 

2p lm ^ Ap 2 



3 2 

(P* + Pim) + ( 2 P™E - 2 Pjm7 )-V + (E 2 + 7 2 ) tP~ = t 2 B 2 - (38) 
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These equations imply that one of the parameters (e.g., E) will be determined once the parameters, 7, e, p r 
and pi m are chosen. 



B. Linear Stability Analysis for the Trimer Case 



We again consider the nonlinear PT-symmetric trimer model with linear and nonlinear gain/loss and examine the 
linear stability of its solutions given in Eq. ([26]) for the solutions given in the previous section. We begin by positing 
the linearization ansatz: 



u(t) = e lEt [a + pe M + Pe x '% v(t) = e im [b + qe M + Q 



(t) = e lht [c + re xt + Re x *] 



(39) 



where p, P, q, Q, r, R are perturbations to the solutions of interest. Plugging Eq. (f39|) into Eq. ([26]) and truncating at 
the linear order in p, P, q, Q, r and R, we derive the following eigenvalue problem: 



AY = iXY 

where Y = (p, — P*, q, — Q*, r, — R*) T and A is the (6 x 6) matrix: 
/ 



an 


-a 2 (p r - ipi m ) 


— e 








[pr + iPim) 


a-22 





e 





— € 





E-2\b\ 2 
-(b*f 


b 2 


— e 





e 


-E + 2\b\ 2 











—e 





055 











e 









e 

-C 2 (p r + ip im ) 
066 



where 



022 



an = E + £7 + 2\a\ 2 (p r - ip im ) 
-E +£7 - 2|a| 2 (p r + ipi™) 
a 55 = E + 2\c\ 2 (p r + ipim), 
E - £7- 2|c| 2 ( ( o r - £pi m ). 



^66 



(40) 



(41) 



(42) 



The solution of this 6x6 eigenvalue problem (and whether the corresponding eigenvalues A possess a positive real 
part) will determine the spectral stability properties of the solutions of the nonlinear-PT-symmetric trimer. 



C. Numerical Results for the Trimer Case 

• Trimer Case I: 

The numerical results for the symmetric solutions of the nonlinear-PT-symmetric trimer (Case I) are shown in 
Fig. 21 Fig. [5] and Fig- [HI with similar notations as in the dimer case. Solutions are found by numerically solving 
Eq. (|34l) . A typical example of the branches that may arise in case I of the trimer is shown for the parameters 
e = 1, E = 1, p r = — 1 and pi m = 1. In this case, we find three branches in the considered interval of parameter 
values. There are two branches which exist up to the point 7 = 2.59 where they collide in a saddle-node 
bifurcation. One of these, the red diamonds branch, is mostly unstable except for 7 € [1-26, 1.33] U [2, 2.11]. For 
7 < 1.26, this branch has one real and one imaginary pair, which become both imaginary for 7 > 1.26 until they 
collide for 7 = 1.33 and yield a complex quartet, which subsequently splits into two imaginary pairs for 7 = 2 
and finally into one real and one imaginary pair for 7 > 2.11. The other one, the black squares branch, is always 
unstable due to one real and one imaginary pair. When these two modes collide, a collision arises between both 
their real and their imaginary (respective) eigenvalue pairs. 
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Aside from the other two branches, the branch associated with the blue stars emerges from 7 = 1 and persists 
beyond the above critical point (and for all values of 7 that we have monitored). In our case, this branch is only 
stable for 7 < 1.25, at which two pairs of imaginary eigenvalues collide and lead to a complex quartet, which 
renders the branch unstable thereafter. This branch behaves very similarly as the one in the linear PT trimer 
case reported in [13j . Both of them bifurcate from zero amplitude after a certain value of 7, persist beyond the 
linear PT critical point and have similar stability properties. 
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FIG. 4: The symmetric solution profiles of Case I in the nonlinear-PT-symmetric trimer with e — 1, E = 1, p r = —1 and 
Pim = 1. The three branches are denoted by blue stars, red diamonds and black squares and their amplitudes (top left), phases 
(top right), real part (bottom left) and imaginary part (bottom right) of the corresponding eigenvalues are shown. See also the 
relevant discussion in the text. 
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FIG. 5: The spectral plane of the linear stability analysis for the symmetric solutions of Case I with e = 1, E = 1, p r = — 1 
and pim = 1, for three different values of 7 = 0.5, 1.5 and 2.5. Each branch is associated with three eigenvalue pairs one of 
which is at due to symmetry. 



To monitor the dynamical evolution of the different branches, we used direct numerical simulations illustrated 
in Fig. [6] for the case of 7 = 1.5. Two of the branches, the blue stars of the left panel and the red diamonds of 
the middle one are oscillatorily unstable for this value of 7, while the black squares branch is always unstable 
due to a real eigenvalue pair. The latter has been found to generically cause the unbounded gain of at least 
one node within the trimer. The oscillatory instability, on the other hand, in the case of the blue branch and 
for 7 = 1.5 can be observed to lead to a longdived periodic exchange of "power" between the three sites. On 
the other hand, for the red diamonds branch of the middle panel, while there is an intermediate stage of power 
oscillations between the three nodes, the ultimate fate of the configuration favors the unbounded growth of at 
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least one node (in fact, two nodes in the example shown) of the trimer. 




t t t 

(a) blue stars branch (b) red diamonds branch (c) black squares branch 

FIG. 6: The time evolution plots of the trimer case I with e = 1, E = 1, p r = —1 and pi m — 1 when 7 = 1.5. For each branch, 
the blue line denotes the nonlinear loss/linear gain site, the red line denotes the nonlinear gain/linear loss site, while the green 
represents the inert site between the two. Notice the oscillatory evolution of the blue solid branch, while the red diamonds and 
black squares lead to ultimate unbounded increase of at least one site within the trimer. 



• Trimer Case II: 

According to Eq. (j35|) . one of the parameters should be determined by the others. We hereby set e = 1, p r = — 1 
and pi m = 1, then E is obtained self-consistenly for a given choice of 7. The solution profiles, which are obtained 
by solving Eq. (|35p numerically, are plotted in Fig. [7J There are three pairs of branches, i.e., six branches of 
solutions in total found in this case. Only one out of each pair is shown in Fig. [TIGS (to avoid cluttering of the 
relevant figures), namely the blue stars, red diamonds and black squares. The other three branches are mirror 
symmetric to these branches, respectively. For example, the existence profile of the mirror symmetric branch 
of the blue stars would be identical to the blue stars shown in Fig. while its stability plot would be mirror 
symmetric about the imaginary axis to the blue stars shown in Fig. [5] Among the three branches shown in 
Fig. [3 two of them emerge at 7 = 2 and persist throughout the range of 7 values considered. It should be 
noticed that the amplitudes are different within each branch. The blue stars (dynamically stable) branch has 
a large A and small B and C, while the red diamonds and black squares branches have a fairly small A and 
large values of B and C . In fact, precisely at the critical point of the branches' emergence, the blue stars and 
the black squares are exact mirror images of each other (i.e., they have the same amplitude for B and the one's 
A is the other's C -and vice versa-). This mirror symmetry is in fact directly reflected in the eigenvalues of the 
linearization around the two configurations, one set of which (for the blue stars) possesses negative real parts, 
while the other (black squares) has mirror symmetric positive ones. As can be perhaps intuitively anticipated, 
the more stable configuration is the one having large amplitude at the nonlinear loss/linear gain site. The third 
branch (red diamonds) is also highly unstable and emerges out of a bifurcation at 7 = 2.05 (to which we will 
return when discussing case III). 

Fig. [5] shows the eigenvalues of the three branches clearly illustrating the fact that they are not symmetric about 
the imaginary axis. This can once again be justified by the asymmetry of the configurations of case II which, 
in turn, break the PT symmetry of the linearization matrix and hence lead to asymmetric spectra. The blue 
stars branch is always stable, as indicated above, and the other two branches are always unstable as 7 increases. 
For instance, for the red diamond branch, there exists (in addition to a zero eigenvalue) an imaginary pair, a 
complex conjugate pair (with a positive real part) and a real eigenvalue. In the case of the black squares branch, 
there are (in addition to the zero eigenvalue) three real pairs (two positive and one negative) and a complex 
conjugate pair (with positivre real part). 

Fig. [5] shows the dynamical plots of the three distinct branches of solutions. The blue stars branch clearly 
preserves its configuration due to its dynamical stability, while for the two unstable branches, their evolution 
gives rise to asymmetric dynamics favoring the loss of the power in a single site (the nonlinear gain/linear loss 
one), and quickly absorbed by a stable state. The latter state appears to be the mirror symmetric of the red 
diamonds asymmetric branch in both cases. This is indeed also a stable dynamical state of the original stationary 
system of equations. Furthermore, this state is expected to exist based on the symmetry breaking bifurcation 
that we will discuss below as giving rise to the red diamonds branch. As remarked above, the stability properties 
of the mirror symmetric branches are mirror symmetric to the ones shown in Fig. [5) This implies that only the 
mirror symmetric branch of the red diamonds is stable, which is, in turn, consonant with our observation that 
it is a potential attractor for the dynamics for 7 = 3 shown in Fig. [SJ 
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FIG. 7: The solution profile of the trimer case II with e = 1, p r = —1 and pi m = 1. The three branches are denoted by blue 
stars, red diamonds and black squares. These branches start at 7 = 2 (except for the red diamonds branch that is initiated at 
7 = 2.05) and exist even when 7 is large. 
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FIG. 8: The spectral stability plots of the trimer case II with e = 1, p r = — 1 and pi m = 1, illustrating the stability of the blue 
stars' branch and the instability of the other two. 
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FIG. 9: The time evolution plots of the trimer case II with e = 1, p T = —1 and pi m = 1 when 7 = 3. The two unstable branches 
(red diamonds and black squares) tend to a dynamically stable configuration which is a mirror image of the red diamonds 
branch. 
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• Trimer Case III: 

Finally, we turn to a consideration, using the same parametric setting as in case II, of the numerical results 
by solving Eq. (|3"5|) for case III in Figs. [TUHl^l Four distinct branches of solutions are observed in this case. 
The branches denoted by red diamonds and black squares exist only for small values of the linear gain/loss 
parameter 7, are stable and terminate at 7 = 0.65. The other two branches, namely the blue stars and the green 
circles collide and terminate at 7 = 2.1. The green circles branch is unstable in this case, due to a complex 
quartet of eigenvalues (observed in Fig. ITTj). On the other hand, the blue stars' branch is stable up to 7 = 2.05 a 
critical point at which branches of case II (the red diamonds branch referred to in case II as having a pitchfork 
bifurcation at the same value and its mirror symmetric image) emerge. Notice that this detail is not discernible 
in the eigenvalue plots of Fig. [TO] 
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FIG. 10: The solution profile for the nonlinear-PT-symmetric trimer of Case III (special symmetric solutions) with e = 1, 
p r = —1 and pi m = 1. The norms are not squared in the top left panel to improve the visibility of the branches (given the 
disparity of the relevant amplitudes). The four branches are denoted by blue stars, red diamonds, black squares and green 
circles. The blue stars, red diamonds and black squares branches always have two pairs of purely imaginary eigenvalues, while 
the green circles branch always has a complex quartet. The blue stars branch terminates with the green circles at 7 = 2.1, 
while the red diamonds and black squares terminate together at 7 = 0.65. 
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FIG. 11: The plots of the spectral plane of the linear stability eigenvalues for the nonlinear-PT-symmetric trimer Case III 
(special symmetric solutions) with e = 1, p r = — 1 and pi m = 1. 



The dynamics of the different configurations are shown in Fig. [12] The blue stars, red diamonds and black 
squares special symmetric branches of solutions are stable and thus preserve their shape. On the other hand, 
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the green circles for 7 = 0.5 are subject to the oscillatory instability predicted by the linear stability analysis. 
This, in turn, results into long-lived oscillatory dynamics of the system, as indicated in the bottom right of 

Fig.ru 
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FIG. 12: The time evolution plots of the nonlinear-PT-symmetric trimer in Case III of special symmetric solutions with e = 1, 
p r — — 1 and pi m — 1 when 7 = 0.5. The only unstable configuration is the green circle one of the bottom right which leads to 
long-lived oscillatory dynamics. 



IV. THE SPECIAL CASE OF LINEAR-PT-SYMMETRIC OLIGOMERS 



As a special case with p r — — 1 and pi rn = of what we did above, the linear PT-symmetric oligomer case example 
has been addressed in the earlier work of [l3j . 

Fig. 1131 shows the profile of the two branches, which are analogues of the nonlinear-PT-symmetric dimer displayed 
in Fig. [TJ The branch denoted by dashed line corresponds to the blue stars branch of case I within the nonlinear- 
PT-symmetric dimer with (— ) sign in Eq. (|11[) and is stable when 7 2 < k 2 — E 2 /4, whereas the solid line branch 
corresponding to the red diamonds branch (of case I in the nonlinear-PT-symmetric dimer) is always stable. The 
linearization around these branches can be performed explicitly in this simpler linear-PT-symmetric case yielding the 

nonzero eigenvalue pairs ±2i\j2(( 2 — 7 2 ) — E^Jl 2 — 7 2 for the first and ±2i^2(e 2 — j 2 ) + E^/t 2 — j 2 for the second 

(notice that the latter can never become real). 

It is relevant to note here that the two branches "die" in a saddle-center bifurcation at 7 = e, as shown in the figure. 
This point coincides with the linear PT-symmetric dimer phase transition. As indicated before, in the nonlinear dimer, 
the two branches die when the restriction (|12l) is no longer satisfied. Nevertheless, the nonlinear solutions of the latter 
case can generally exist past the linear phase transition (and even arbitrarily past that as in the case II solutions) 
and moreover asymmetric solutions can exist due to the interplay of linear and nonlinear gain/loss a feature absent 
in the simpler linear-PT-symmetric dimer. 

As an analogue of the case I solutions of the nonlinear-PT-symmetric trimer, we present a prototypical example of 
the branches that may arise in the case of the linear-PT-symmetric trimer in Fig. fT?]for E = e = 1. There are three 
distinct branches. Two of them collide in a saddle-center bifurcation (for 7 = 1.043) and disappear thereafter. The 
other one emerges and bifurcates from zero amplitude for 7 > \/2e 2 — E 2 and persists beyond the critical point of 
the linear PT phase transition jpx = \/2e, (although it is unstable in that regime). Hence the feature of solutions 
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FIG. 13: (Color online) The two branches of solutions for the dimer problem are shown for parameter values e — E — 1. (a) 
The amplitude of the sites, (b) their relative phase, and (c) the (nontrivial) squared eigenvalue of the two branches. The solid 
line corresponds to the alw ays stable branch, while the dashed line corresponds to the branch, which acquires a real eigenvalue 
pair above a certain 7 = ^/ e' 2 — E 2 /4. Reprinted with permission from 
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FIG. 14: (Color online) Existence and stability of solutions for the case of the trimer analogous to Fig. [13] with parameters 
E = 1 and 4>t normalized to (without loss of generality) . There are three branches: uP^ (dash-dotted line), m' 2 ' (dashed line), 
and it' 3 ' (solid line). For each branch, two curves in (a) stand for A,B,C (since C = A), and two curves in (b) stand for <f) a 
and 4> c . Panels (c) and (d) are the illustration of both real and imaginary parts for the eigenvalues. There exist two solutions 
(the dashed and the dash-dotted) one of which is stable (dashed) and one generically unstable (dashed-dotted) which collide 
and disappear for 7 = 1.043. The third solid branch emerges from the limit 7 = 1 and persists thereafter (for all values of 7 
considered), although it becomes unstable for 7 > 1.13. Reprinted with permission from [l3T |. 
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persisting past the linear PT phase transition exists even in the linear-PT-symmetric trimer, yet other more complex 
features do not appear in this setting. A canonical example thereof is the existence of asymmetric solutions (which, 
in turn, possess asymmetric spectra). The latter trait is amply present in the nonlinear-PT-symmetric trimer of the 
previous section. 

It should be added that in addition to the oligomers of the dimer and the trimer variety, recently there has also been 
considerable interest towards the study of quadrimcr settings. Such settings were again initiated in the examination 
of [IH, where the particular (and simpler) case of (—47, —17, +17, +17) i.e., a linear-PT-symmetric quadrimer with 
two lossy sites on the one side and two gain nodes on the other side was considered. The considerations of [l3[ were 
generalized in the very recent study of [26J, which examined the configuration (—171, —172, +«72, +*7i) he-, the most 
general bi-parametric gain/loss family of quadrimers. In both cases, the possibility of solutions that emerge purely 
in the nonlinear regime and of ones that continue past the critical point of the underlying linear PT phase transition 
was identified. Another reason why such quadrimers are of interest is that they can be thought of as the prototypical 
building blocks (if placed on a square "plaquette" ) of a two-dimensional PT-symmetric lattice [30] . To the best of our 
knowledge, nonlinear-PT-symmetric quadrimers, along the lines considered herein have not been examined to date. 



V. CONCLUSIONS AND FUTURE CHALLENGES 



In the above study, we illustrated some interesting characteristics which emanate from the interplay of nonlinearity 
with PT-symmetric linear Hamiltonians in the case of oligomer configurations. The basic underlying premise which 
has been explained in the recent works of [27], H(| and has been explored in numerous others such as [13|, H3, H3, H3] 
is that the nonlinear and the PT-symmetric-linear part of the system at hand no longer commute and hence give 
rise to novel phenomenology that is not expected for linear PT symmetric systems. A principal element of this 
phenomenology that arises even for PT-symmetric oligomers (trimers, quadrimers, etc.) with merely linear gain/loss 
is the fact that nonlinear solutions may exist that do not have a corresponding linear limit and which, in fact, defy 
the threshold for the linear PT phase transition in that they exist for arbitrary gain/loss parameter values past 
that critical point. The introduction of a nonlinear gain/loss pattern considered herein for dimers and trimers (and 
earlier in a more cursory way in [24j for dimers) presents additional possibilities stemming from the interplay of linear 
and nonlinear gain/loss profiles. These include among others the emergence of asymmetric solutions which not only 
are involved in symmetry breaking (pitchfork) bifurcations but also produce asymmetric linearization matrices with 
spectral properties that reflect this asymmetry. 

We believe that this direction of studies is particularly intriguing for further progress, especially as the complexity 
of the problem increases within the confines of a full one-dimensional chain, but also even for fundamental two- 
dimensional entities, such as the quadrimer based plaquettes. These themes constitute pristine territory for further 
exploration at the theoretical level. Naturally, on the other hand, a potential generalization of the earlier experiments 
of |5j towards the inclusion of nonlinear amplification and amplitude-dependent loss in balance with each other would 
be most worthwhile to consider in order to take advantage of the considerable additional wealth of phenomenology of 
the latter system. 
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